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Abstract 

We report a preliminary study of the branching fractions and cp' distributions of exclusive charm- 
less semileptonic B decays, using events tagged by fully reconstructing one of the B mesons in a 
hadronic decay mode. These results are obtained from a data sample that contains 535 x 10® BB 
pairs, collected near the T(45) resonance with the Belle detector at the KEKB asymmetric energy 
e^e" collider. 

PACS numbers: 13.25.Hw, 13.30.Ce, 14.40.Nd 
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The Standard Model (SM) of particle physics contains a number of parameters whose 
values are not predicted by theory and must therefore be measured by experiment. In 
the quark sector, the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix govern 
the weak transitions between quark flavours, and precision measurements of their values 
are desirable. In particular, much experimental and theoretical effort is currently being 
employed to test the consistency of the CKM formalism 0], by examining the Unitarity 
Triangle most relevant to the decays of B mesons. 

The precision to which the angle sin20i characterising indirect CP violation in 6 — > ccs 
transitions has improved to approximately 5% 0. This makes a precision measurement 
of the length of the side of the Unitarity triangle opposite to (pi particularly important 
as a consistency check of the SM picture. The length of this side is determined to good 
approximation by the ratio of the magnitudes of two CKM matrix elements, |Ktb|/|Ki:fe|- 
Both of these can be measured using exclusive semileptonic B meson decays. Using charmed 
semileptonic decays, the precision to which \Vcb\ has been determined is of order 2%. On 
the other hand \Vub\, which can be measured using charmless semileptonic decays, is the 
most poorly known of the CKM matrix elements. Both inclusive and exclusive methods of 
measuring \ Vub\ have been pursued, with the inclusive methods giving a value to a precision 
of 7-8%. The exclusive determination of \Vub\ currently has a precision poorer than 10%. 
It is the aim of an ongoing programme of measurements at the B factories to improve this 
precision to better than 5%, for comparison with the inclusive results, which have somewhat 
different experimental and theoretical systematics. This would provide a sharp consistency 
test with the value of sin20i. 

In this paper we present preliminary studies of the exclusive semileptonic decays B — > 
TT+^z/, B TT^iu, B P^C-V) B — >■ p^^v and B uoCv using a full reconstruction tagging 
technique to identify candidate B mesons. This measurement is based on a data sample that 
contains 535 xlO^-Bi? pairs, collected with the Belle detector at the KEKB asymmetric- 
energy e"*"e^ (3.5 on 8 GeV) collider KEKB operates at the T(4S') resonance (i/i = 
10.58 GeV) with a peak luminosity that exceeds 1.6 x 10^^ cm~^s~^. The T(4S') is produced 
with a Lorentz boost of /57 = 0.425 nearly along the electron beamline {z). 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
sample of 152 xlO^ BB pairs, while a 1.5 cm beampipe, a 4-layer silicon detector and a 
small-cell inner drift chamber were used to record the remaining 383 xlO^-B-B pairs|5,]- 

Table HI lists all current measurements of branching fractions for exclusive B X^Ci' de- 
cays, where denotes a light meson containing a u quark. This table is based on a recent 
compilation by the Heavy Flavour Averaging Group (HFAG) 0. Three methods of iden- 
tifying signal candidates have been employed in these measurements, denoted "untagged", 
"semileptonic tagged" or "full reconstruction tagged" . Only BaBar have reported results to 
date based on the third of these l^l 13 • 



The most precise measurements at the present time come from the untagged analyses of 
CLEO [7j and BaBar As more integrated luminosity is accumulated by the 5-factory 
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experiments, the full reconstruction tagging technique will become the most precise method. 
It holds the advantage of providing the best signal to background ratio, offset by the lowest 
efficiencies. 

TABLE I: Measurements of branching fractions of exclusive B Xu^u decay modes. In each 
case, the first error is statistical, the second experimental systematic, the third due to form factor 
uncertainties for the signal mode, and the fourth, when present, due to form factor uncertainties 
for crossfeed modes. U indicates untagged method, S semileptonic tagging method, and F full 
reconstruction tagging. 



Experiment 




Mode 


BB 


Branching Fraction 


and Tag Method 






r-i r\fi1 


r-i rv — 41 

[10 


CLEO [7] 


U 


5° - 




9.7 


1.33 ± 0.18 ± 0.11 ± 0.01 ± 0.07 


BaBar [8] 


U 


- 




86 


1.38 ± 0.10 ± 0.16 ± 0.08 


Belle [9] 


S 


B^ - 




275 


1.38 ± 0.19 ± 0.14 ± 0.03 


Belle [2j 


S 


B+ - 




275 


0.77 ± 0.14 ± 0.08 ± 0.00 


BaBar [10] 


s 


B — 


— ■ /) 

^ IT tU 


232 


1.03 ± 0.25 ± 0.13 


BaBar [11] 


s 


B+ - 


ir^iu 


88 


1 on _|_ n QV _|_ n OO 


BaBar [12] 


F 


B^- 




233 


1.14 ±0.27 ±0.17 


BaBar [12] 


F 


B+ - 




233 


0.86 ±0.22 ±0.11 


BaBar [12] 


F 


B - 




233 


1.28 ±0.23 ±0.16 


CLEO [T] 


U 


B^- 


^ p-lv 


9.7 


2.17 ±0.34l[}:^J ±0.41 ±0.01 


CLEO [13] 


U 


B^- 


^ p-lv 


3.3 


2.69 ± 0.4ll|j ± 0.50 


BaBar [14] 


F 


B^- 


^ p-lv 


88 


2.57 ±0.52 ±0.59 


BaBar [15] 


U 


5°- 


> p- ev 


55 


3.29 ±0.42 ±0.47 ±0.60 


BaBar [8] 


U 


B^- 


^ p-lv 


83 


2.14 ±0.21 ±0.51 ±0.28 


Belle [2] 


S 


B^- 


^ p-lv 


275 


2.17 ±0.54 ±0.31 ±0.08 


Belle [9] 


S 


B+ - 


p^lv 


275 


1.33 ±0.23 ±0.17 ±0.05 


Belle [16] 


U 


B+ - 


ujIv 


85 


1.3 ±0.4 ±0.2 ±0.3 


CLEO [7] 


U 


B+ - 


r]£v 


9.7 


0.84 ±0.31 ±0.16 ±0.09 



In this analysis we fully reconstruct one of the two B mesons from the T(45') decay (-Btag) 
in one of the following hadronic decay modes, B~ D^*^^7i~, B~ D^*^^p~, B~ — > D^*^^ai, 
B- B° DW+tt", 5° ^ D^*^+p-, 5° D(*)+ar or 5° ^ !)(*)+/}(*)- 

[l7l |. Decays are identified on the basis of the proximity of the beam-energy constrained 
mass Mbc = y'(-^bcam)^ — {pbY and AE = Eb — -Ebeam to their nominal values of the 
B meson rest mass and zero, respectively. Here -Ebeam, Pb and Eb are the beam energy 
and the measured 3-momentum and energy of the Stag candidate in the T{4S) rest frame 
respectively. If multiple tag candidates are found, the one with values of Mbc and AE closest 
to nominal is chosen. Events with a i?tag satisfying the selections Mbc > 5.27 GeV/c^ and 
—0.08 < AE < 0.06 GeV are retained. The charge of the -Btag candidate is necessarily 
restricted to Qtag = or Qtag = ±1 by demanding that it is consistent with one of the above 
decay modes. 

Reconstructed charged tracks and ECL clusters which are not associated with the i?tag 
candidate are used to search for the signal B meson decays of interest recoiling against the 
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Stag- Photons identified with isolated ECL clusters which have a laboratory energy of less 
than 50 MeV are ignored. 

Electrons are identified using information on dE/dx from the CDC, response of the ACC, 
shower shape in the ECL, and the ratio of the energy deposited in the ECL to the momentum 
determined from tracking. The signals in the KLM are used to identify muons. Charged 
kaons are identified based on the dE/dx information from the CDC, the Cerenkov light yields 
in the ACC and time-of-fiight information from the TOF counters. Any charged particles 
which are not identified as leptons or kaons are taken to be pions. 

Photons whose direction in the laboratory frame lies within a 5° cone of the direction of 
an identified lepton are considered to be bremsstrahlung. The 4-momentum of the photon 
is added to that of the lepton and the photon is not considered further. 

Neutral pions are reconstructed from pairs of photons whose invariant mass lies in the 
range [0.120, 0.150] GeV/c^, of order ±3(7 of the 7r° mass. Charged p meson candidates are 
reconstructed via the decay p'^ tt"^ tt^ where the invariant mass of the pair of pions is 
required to lie in the range [0.570,0.970] GeV/c^. Neutral p meson candidates are similarly 
reconstructed from pairs of oppositely charged pions, with the requirement that m7r+7r- is 
in the range [0.703,0.863] GeV/c^. Finally, u; candidates are reconstructed from a; — > 7r+ 
7r~ 7r° with m^+^-7rO in the range [0.690,0.850] GeV/c^. In events where more than one 
hadron candidate, denoted X„, of a given type is identified amongst the recoil particles, the 
candidate with the highest momentum in the T{4S) rest frame is chosen. 

To isolate signal candidates, several requirements are placed on the recoil system. There 
must be one lepton candidate present only, with pi^p > 0.4 GeV/c in the T(45') rest frame. 
The total charge of the recoil system, Qrccoib is required to be if a neutral tag has been 
identified, and ±1 if a charged tag has been identified. In the charged case, the sign of Qrccoii 
must be opposite to that of Qtag- In the neutral case, we do not make any requirement on 
the sign of the lepton charge with respect to the Stag, to allow for mixing. 

The number of charged recoil particles is required to correspond to one of the sought 
modes, i.e. one for B Ti^iu (the lepton), two for B vr+fz/ and B — > p'^lv (the lepton 
plus one charged pion) and three for B — > p^lv and B uj£u (the lepton plus two charged 
pions). Additionally, the number of recoil 7r° candidates is required to be consistent with 
one of the sought modes. In order to increase efficiency, however, we allow more than the 
necessary number in some cases: we require no vr'' candidates to be present for B — >■ Tr'^ii' 
and B — > p^iu modes, and at least one 7r° for the B n^iu, B p'^iu and B uoiu 
modes. Additionally, we require that there be no more than 0.3 GeV of residual neutral 
energy present on the recoil side in the T(45') rest frame, after any photons contributing to 
the Xu candidate have been removed. 

Signal events are identified by examining the missing mass squared (M^j^g) distributions. 
If the tagging B is correctly reconstructed and the correct lepton and hadron candidate 
have been identified on the recoil side, then (ideally) all missing 4-momentum is due to the 
remaining unidentified neutrino. The square of the missing 4-momentum for signal events 
should therefore be close to zero, and applying this requirement provides a very strong 
discrimination between signal and background. 

Background contributions come from several sources. These include semileptonic decays 
resulting from h clu transitions, denoted B — > Xciv, which have significantly larger 
branching fractions to the channels under study; continuum + qq processes; and 

cross feed from one B X^^f channel into another. The contributions of these backgrounds 
are studied using Monte Carlo (MC) simulated data samples generated with the EvtGeb 
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package [13]. Generic BB and continuum MC samples equivalent to approximately twice 
the integrated luminosity of the real data set are used. The model adopted for B — ^ D*ii' 
and B Diu decays is based on HQET and parametrisation of the form factors , while 
B D**iiy decays are based on the 1SGW2 modeljii. A non-resonant B —>■ D^*^nii' 



component based on the Goity-Roberts prescription 2^ is also included 

A separate MC sample equivalent to approximately seven times the integrated luminosity 
of the real data set is used to simulate the signal channels and crossfeed from other B — > X„£z/ 
decays. Models for the exclusive modes are based on Light Cone Sum Rules (LSCR) for tt 
l^l, p and u ji^ modes and ISGW2 ji^l for other exclusive modes. 



Radiative effects associated with the lepton and resulting from high er-order QED pro- 
cesses are modelled in all MC samples using the PHOTOS package 124\. All generated MC 
events are passed through a full simulation of Belle detector effects based on GEANT 3.21 
0. 

Figure □ shows the observed M^jgg distributions for the five decay modes. Also shown in 
this figure are the fitted signal component, the fitted uiiy crossfeed, and the fitted contri- 
bution from other backgrounds, which is dominated by S — >■ Xf.ii' decays. The shapes of 
these components are taken from MC and the normalizations are fit parameters. The fitting 
method follows that of Barlow and Beeston |2^ and takes into account finite Monte Carlo 
statistics. The fitted event yields obtained are 48 ± 8 for the B —>■ Tr^iu mode, 35 ± 7 for 
B Ti'^iu, 41 ± 9 for 5 ^ p+iu, 61 ± 9 for 5 p'^iu and 27 ± 9 for 5 ^ uiu. 

We extract the partial branching fractions in bins of g^, the invariant mass squared of the 
lepton-neutrino system, in order to minimise the systematic error which arises from the lack 
of precise knowledge of the shape of the form factors. These are shown in Figure„ Three bins 
of are chosen, commensurate with available statistics, to 8 GeV^/c^, 8 to 16 GeV^/c^, 
and greaterthan 16 GeV^/c^. The neutrino 4-momentum is determined from the missing 4- 
momentum vector using p^, = (IPmissl^Pmiss), where Pmiss is the missing 3-momentum vector 
in the T(4S') rest frame. The resolution obtained varies from 0.20 GeVVc^ for the 
B TT^iu channel to 0.27 GeV^/c^ for the B coiu channel. 

Table |n] summarises the result of a preliminary study of the contributions to the total 
systematic error for the branching fractions summed over the three bins, for the B — >■ 
TT and B ir^iu modes. These are broken down into the following categories; those 
arising from detector simulation, such as charged track reconstruction efficiency, particle 
identification and neutral cluster reconstruction; uncertainties in the luminosity; and effects 
of the form factor models used and assumed branching fractions in the MC. 

The effects of model dependence of the form factor shapes assumed in the B —>■ X^^v MC 
used for signal efficiency and crossfeed background estimates have been studied by comparing 
the fitted yields obtained using the default model implemented in the MC, which is LCSR 
(23I I 2^ 1 ■ and the ISGW2 model ji^. This is achieved by reweighting the MC events on 
an event-by-event basis based on their generated values of and angular variables. The 
variation between these two models in predicting the shapes of the distributions for the 
pseudoscalar and vector modes typifies the spread between available models for the dynamics 
of these decays. 

Effects due to the uncertainties in the branching fraction normalisations for h —>■ uiu and 
b — s> ciu decays in the signal and background MC samples were studied by varying in turn 
the B n+iu, B n^iu, B p+iu, B p°£z/, B ^ uiu, B D+iu, B D'^iu, 
B — s> D*^£i' and B D*^ii> branching fractions by their measurement errors as quoted 
by the Particle Data Group ji^. A reweighting technique is again used, and fitted yields 
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with and without reweighting are compared. The maximum observed spread in the fitted 
branching fraction is assigned as systematic error. 

The effects of finite MC statistics are taken into account in the fitting procedure [i^ 
and are reflected in the errors on the obtained branching fractions. Since the available MC 
samples are rather limited in statistics, variations of the assumptions on form factor shapes 
and normalizations can be absorbed by the present fits to a significant extent. 



TABLE II: Results of a preliminary study of sources of systematic uncertainty. 



Source of error 


Assigned systematic error 




B TT+iu B -K^iu 



Detector Simulation: 






Pion track finding efF. 


1.3% 




vr'^ reconstruction eff. 




4.6% 


Lepton track finding efF. 


1% 


1% 


Lepton identification 


2.1% 


2.1% 


Charged kaon identification 


2.0% 




Combined 


3.3% 


5.2% 


N{BB ) uncertainty 




1.3% 


Form Factor Shapes: 






vr (LCSR ^ ISGW2) 


2.2% 


1.4% 


p,io (LCSR ^ ISGW2) 


0.1% 


2.4% 


Branching Fractions: 






b — > uiu 


0.0% 


2.2% 


b^clu 


0.7% 


2.2% 


Total systematic error 


4.2% 


6.8% 



The partial branching fractions in bins of are given in Table UTTl for the B —>■ Tc^iu and 
B — > n^iu modes, for which reliable preliminary systematic errors have been estimated. The 
systematic errors are included for each bin as well as the sum over bins. Figure |21 presents 
the shapes of the partial branching fractions for all five modes as a function of g^, where the 
statistical and systematic errors have been added in quadrature. The systematic errors for 
the vector modes in this figure should be considered as very preliminary. 



TABLE III: Partial branching fractions in three bins of q'^. These are summed to give the full 
branching fraction quoted in the "Sum" column. Errors are statistical and systematic. 







AB [10-4 


B [10-^ 




Mode 


< g2 < 8 


8 < < 16 




> 16 


Sum 






(GeVVc^) 


(GeVVc^) 




(GeVVc2) 


(GeVVc^) 


B TT+lu 


0.50 ±0.14 ±0.02 


0.68 ±0.18 ±0.03 


0.31 ±0.12 ±0.01 


1.49 ±0.26 ±0.06 


B TT^iu 


0.28 ±0.09 ±0.02 


0.22 ±0.08 ±0.04 


0.36 ±0.12 ±0.02 


0.86 ±0.17 ±0.06 



In summary, we have made a preliminary study of the partial branching fractions as a 
function of for five semileptonic decay channels of B mesons to charmless final states. 
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using a full reconstruction tag method. Summed over the three bins we obtain the 
following estimates of the branching fractions for the pion modes: B{B ^ vr+^i/) = 
(1.49 ± 0.26 ± 0.06) X 10-^ B{B ^ Ti^iu) = (0.86 ± 0.17 ± 0.06) x 10~^ where the first 
error is statistical and the second a preliminary estimate of the systematic error. Whilst the 
statistical precision of these measurements is limited at present, the potential power of the 
full reconstruction tagging method, when it can be used with larger accumulated B-factory 
data samples in the future, can clearly be seen. 
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FIG. 1: Missing mass squared (M^j^g) distributions after all cuts, for (a) B ■K'^iu, (b) B ir^iu, 
(c) B p'^^i^i (d) B p^tu, and (e) B — > ojIv modes. Data is indicated by the points with error 
bars. The blue histogram (lightest shade in greyscale) shows the fitted prediction based on the 
LCSR model [i^l H^- The green histogram (middle shade in greyscale) shows the fitted h — > u^v 
background contribution. The crimson histogram (darkest shade in greyscale) shows the fitted 
background contribution from other sources. The fitting method is explained in the text. 
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FIG. 2: Partial branching fractions as a function of for the five signal modes (a) B — > ir^iv, 
(b) 5 — 7r°^z^, (c) B p'^if, (d) B p^lv, and (e) B — a;£i^. Errors shown are statistical and 
preliminary systematic, added in quadrature. 
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